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We discuss the properties of ghost and gluon propagators in the deep infrared momentum region 
of Landau gauge Yang-Mills theory. Within the framework of Dyson-Schwinger equations and the 
functional renormalization group we demonstrate that it is only a matter of infrared boundary 
conditions whether infrared scaling or decoupling occurs. We argue that the second possibility is 
at odds with global BRST symmetry in the confining phase. For this purpose we improve upon 
existing truncation schemes in particular with respect to transversality and renormalization. 
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I. INTRODUCTION 

In the past decade much progress has been made in 
the understanding of the infrared sector of QCD de- 
scribed in terms of Green's functions. Various meth- 
ods have been employed for this purpose, in particu- 
lar Dyson-Schwinger equations (DSEs) 
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stochastic quantization methods [35|, l36|, effective the- 
ory ap p roaches l37i . 38, 39. l40| and lattice gauge theory 
m m, S, H Si H EtI, m m MMM- However, 
as the full Green's functions are non-perturbative, their 
determination is notoriously complicated, and each and 
every method faces individual technical problems. 

Already in pure Yang-Mills theory the determination 
of these Green's functions poses a considerable challenge. 
It is complicated by the fact that these are in general 
gauge-dependent quantities. However, this fact can also 
be used to simplify matters by an appropriate choice of 
gauge. Most commonly Landau gauge has been used in 
this respect. 

Over the past ten years substantial progress has been 
made using functional methods in this gauge. Based on 
the pioneering works ref. P, 0, [IHl powerful tools have 
been developed to analyze the whole tower of Dyson- 
Schwinger and functional renormalization-group equa- 
tions in the deep infrared 0, As a result a self- 
consistent scenario has been identified which consists of 
a dressed ghost propagator more divergent in the infrared 
than its tree-level counterpart and an infrared suppressed 
gluon propagator with a finite or even vanishing dressing 
function at zero momentum. This solution agrees well 
and in fact supports the Kugo-Ojima confinement sce- 
nario developed in ref. [s^ . As we shall argue below this 
solution is the only one that has this property. 

In lattice gauge theory substantial progress has been 
made in the past years to overcome intrinsic problems 



related to discretization and volume artefacts as well as 
gauge fixing ambiguities. Recent results on very large 
volumes j46l |47| . |49| produce an infrared finite gluon 
propagator and, more important, an infrared finite ghost 
in disagreement with the continuum results discussed 
above. However, the situation on the lattice is still 
not settled yet. These results are affected from dis- 



cretization artefacts [5J] and effects due to Gribov copies 
psi . m, [HI, [HI], to a significant extent. We will discuss 
these issues below in section IVIII 

Nevertheless, the lattice results have inspired a number 
of recent works in the continuum theory aiming at solu- 
tions of Dyson-Schwinger equations with a finite ghost in 
the infrared [13, [13, HI]- In this work we confirm that 
these solutions are present if the functional equations are 
taken on their own correctly. However, as we will argue in 
section lTl Bl these solutions cannot represent the infrared 
behavior of Yang-Mills theory in the confined phase while 
at the same time maintaining global color symmetry and 
BRST symmetry. 

This work is organized as follows. We will discuss 
the basic setting, the non-perturbative regime of Yang- 
Mills theory, and one set of methods to determine these 
Green's functions, functional equations including bound- 
ary conditions, in section[TTl In particular, in section lll Bl 
we will discuss the roles of global symmetries and iden- 
tify boundary conditions necessary to obtain a BRST- 
symmetric, confining solution. The derivation of the 
functional equations requires a discussion of certain for- 
mal aspects, renormalization and gauge invariance done 
in section [HTI In section ITVl we then discuss the impact of 
the boundary conditions on the infrared behavior of the 
ghost and gluon propagators, the simplest Green's func- 
tions. As an example we reconsider a truncation scheme 
for the DSEs which has been introduced in [3|. For the 
ghost DSE this truncation is identical to those chosen by 
Boucaud et al. [131 and Aguilar et al. We also intro- 
duce a modified scheme which maintains transversality 
in the gluon DSE and is free of quadratic divergencies. 
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We solve the coupled system of ghost and gluon DSEs 
with different choices of boundary conditions and verify 
the appearance of two different types of solutions in the 
infrared. We compare our approach with others and com- 
ment on various technical issues. In section IVTl we repeat 
this exercise in the framework of the functional renor- 
malization group. We compare and comment on recent 
lattice results in section IVlIl We also discuss in sectionllll 
and IVIII the role of different non-perturbative extensions 
of the Landau gauge. In addition, we show in section lVIIII 
also that neither type of solution can be associated with a 
massive gluon characterized by a gauge-independent pole 
mass. Finally, in section IIXI we summarize and discuss 
our results further, the main result being an understand- 
ing of how and why both types of solutions exist, and, in 
the framework of local quantum field theory, which could 
be the preferred one. 

II. NON-PERTURBATIVE YANG-MILLS 
THEORY 

A. Gauge-fixing 

The starting point for the derivation of functional 
equations is the gauge-fixed path-integral. Here we shall 
use Landau gauge, 

d,A; = o. (1) 

Within the standard Faddeev-Popov gauge fixing proce- 
dure this leads to the gauge fixed action in Euclidean 
space-time 

S^J d'x {F^,F^^, + -c'-d.Dfc') , (2) 

with the ghost fields c° and c° and 

= d,Al~d.A;+gr'^AlA: 
Df = 5^'d,,+gr'Al. (3) 

This procedure is only well-defined within perturbation 
theory. In general there are several gauge equivalent con- 
fingurations satisfying the gauge condition ([1]) . These are 
called Gribov copies [53, [S^. Indeed, the related path in- 
tegral vanishes identically, as contributions from different 
gauge copies cancel out[59[. 

The latter problem can be resolved by restricting the 
configuration space to the first Gribov region js^]. This 
region is characterized by the condition that the Faddeev- 
Popov operator Af°^ = —d^D"^ is positive semi-definite. 
This gauge prescription is called minimal Landau gauge. 
However, it is still not sufficient to single out a unique 
gauge copy [6l|. Furthermore it breaks global BRST- 
invariance, since global ERST transformations relate dif- 
ferent Gribov copies [s^. 

A unique selection of gauge copies can be achieved 
by restricting the configuration space to the fundamen- 
tal modular region. This region is characterized by the 



global minimum of the gauge fixing potential |6l| . 

HA)^c j d^xA^AI, (4) 

with c a positive constant. The corresponding gauge is 
also called the absolute Landau gauge [56J . Based on the 
work [3^ one may hope that absolute Landau gauge is 
compatible with global ERST symmetry in the thermo- 
dynamic limit. This gauge is discussed further in section 
IVHl 

B. Global symmetries and confinement 

The appearance of Gribov copies is related to the lo- 
cal (i.e. perturbative) nature of the standard Faddeev- 
Popov procedure without the supplemental condition of 
minimising Consequently, to obtain a well-defined 
non-perturbative functional integral, it is necessary to 
improve upon this procedure. 

A promising approach in this respect is to upgrade 
the gauge fixing procedure to a globally defined (i.e. 
non-perturbative) BRST-symmetry that admits global 
ERST-charges. Then, the path integral is non-vanishing 
and well-defined. This is also necessary for a successful 
implementation of ERST symmetry on the lattice [59| . 
Contrary to the perturbative setting such a procedure 
implies a weighting of Gribov copies such that the con- 
tributions from different gauge copies do not cancel out. 
Very recently within the framework of topological field 
theories a globally valid ERST quantisation has been put 
forward in Ref. [6(J |. 

Functional methods such as Dyson-Schwingcr or func- 
tional Renormalization Group equations have the advan- 
tage that such an explicit construction of a gauge fixing in 
terms of a topological field theory can be avoided. The 
reason is that global constraints from the gauge fixing 
procedure do not alter the form of the functional equa- 
tions but constrain the boundary conditions of Green's 
functions at vanishing momenta in the infrared. The re- 
lated infrared constraint has been deduced by consistency 
arguments for a local covariant quantum field theory, 
see Q. This constraint is directly related to the global 
ERST-charge: the Kugo-Ojima confinement scenario is 
based on a formulation with well-defined global ERST- 
charges [S^l- These enable the definition of the physi- 
cal part Hphys of the state space of Yang-Mills theory 
in terms of the cohomology of the ERST charge oper- 
ator and is the basis of the well-known ERST-quartet 
mechanism. Given furthermore that one is interested 
in a confining solution of infrared Yang-Mills theory a 
well-defined global color charge is then necessary to en- 
sure that Tiphys contains only colorless states. In Landau 
gauge this enforces an infrared enhancement of the ghost 
propagator [53|; the ghost dressing function is then in- 
frared divergent. 

In functional methods this enhancement can be imple- 
mented as an infrared renormalization condition. This 
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condition leads to a unique [101 (scaling) solution of the 
whole tower of functional equations for the one-particle 
irreducible Green's functions of Yang-Mills theory. Hence 
it implicitly defines the unique gauge fixing with well- 
defined global BRST-charges. 

In turn, given confinement, an infrared solution with 
finite ghost at zero momentum (termed 'decoupling' be- 
low) implies broken global gauge and BRST symme- 
tries. Indeed, all known BRST-quantizations that are 
compatible with an infrared finite ghost even break off- 
shell BRST [11,140, [el. The only possibifity for the de- 
coupling solution to coexist with a globally well-defined 
BRST charge is in a Higgs phase, where the breaking 
of global color symmetry implies the existence of super- 
selection sectors. 

Hence, in the continuum we have a unique scenario 
with global BRST-charges and an infrared enhanced 
ghost. These global constraints supplementing Landau 
gauge are most easily implemented within functional 
methods with consistent renormalization conditions fixed 
in the infrared, as will be demonstrated in detail in sec- 
tion [IVl 



C. Functional relations 

We proceed by briefiy discussing the functional rela- 
tions for the Green's functions we use in the present work 
following [loj . 

A quantum field theory or statistical theory can be 
defined uniquely in terms of its renormalized correlation 
functions. They are generated by the effective action F, 
the generating functional of IPI Green's functions. For 
the functional RG equation we also consider the effective 
action in the presence of an additional scale /c, where the 
propagation is modified via fc-depcndent terms 
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identities are projections of the quantum equations of 
motion, whereas flow equations can be read as differen- 
tial DSEs, or DSEs as integrated flows. Written as a 
functional relation for the effective action F^ the DSE 
reads, e.g. (23 . 
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where the operators (j)op are defined as 



(j){x) , 



and 
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(8) 



is the full field dependent propagator for a propagation 
from 01 to ■ The functional derivatives in ([7|) act on 
the corresponding fields and generate one loop and two 
loop diagrams in full propagators. The functional DSE 
([7]) relates IPI vertices, the expansion coefficients of F^, 
to a set of one loop and two loop diagrams with full prop- 
agators and full vertices, but one classical vertex coming 
from the derivatives of Sc\. We emphasize that the DSE 
([7]) only implicitly depends on the regularization via the 
definition of the propagator in ([8|) . 
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are fc-dcpendent regulator functions. 
Within the standard choice k is an infrared cut-off scale, 
and the functions R cut-off the propagation for momenta 
smaller than k. The regularized effective action F^ is 
expanded in gluonic and ghost vertex functions and reads 
schematically 



F, 



E 



}_ p(2n,m) qu j^r. 



m\n\ 



(6) 



in an expansion about vanishing fields (p = {A, C, C). In 
([6[) an integration over momenta and a summation over 
indices is understood. The effective action F^ satisfies 
functional relations such as the quantum equations of mo- 
tion, the Dyson-Schwinger equations (DSEs); symmetry 
relations, the Ward or Slavnov- Taylor identities (STIs); 
as well as functional RG or flow equations (FRGs). All 
these different equations relate to each other, for a de- 
tailed discussion see [l^]. Indeed, the Slavnov- Taylor 



FIG. 1: Functional DSE for the effective action. Filled circles 
denote fully dressed field dependent propagators ((Sjl. Empty 
circles denote fully dressed field dependent vertices, dots de- 
note field dependent bare vertices. 

For Yang-Mills theory a diagrammatic representation 
of the structure of the functional DSE ([7]) is shown in 
Fig. [T] The rhs is given in powers of the ficld-dcpcndcnt 
fully dressed propagator G^^ [(/)], and its derivatives, as 
well as the field dependent bare vertices. The momen- 
tum scaling of Green's functions is directly related to the 
scaling of these building blocks. 

Wetterich's flow equation [2^ for the effective action 
of pure Yang- Mills theory reads, e.g. [24, 28, 29, 3i,^], 



dtTk\ 



i2n) 



{2n) 



GabmP^p) dtR'^'ip), (9) 
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where t = Ink. The flow ([9]) relates the cut-off scale 
derivative of the effective action to one loop diagrams 
with fully dressed field-dependent propagators. We can 
contrast the diagrammatic representation of the DSE in 
Fig. [1] with that of dH), given in Fig. H 




FIG. 2: Functional flow for the effective action. Filled circles 
denote fully dressed field dependent propagators ((8)1. Crosses 
denote the regulator insertion dtR. 

Fig. [5] shows the structure of the functional flow 
The rhs is given by the ficld-dcpcndcnt fully dressed prop- 
agator G,j,^[(j)\ and the regulator insertion dtR. The stan- 
dard use of ([9]) is to take a regulator function R{p'^) which 
tends towards a constant in the infrared and decays suf- 
ficiently fast in the ultraviolet, and hence implements an 
infrared cut-off. 




FIG. 4: Functional RG equations for the gluon and ghost 
propagator. Filled circles denote dressed propagators and 
empty circles denote dressed vertex functions. Crosses de- 
note insertions of the regulators. 





FIG. 3: Dyson-Schwinger equations for the gluon and ghost 
propagator. Filled circles denote dressed propagators and 
empty circles denote dressed vertex functions. 

The diagrammatical expressions for the exact Dyson- 
Schwinger equations (DSEs) for the ghost and gluon 
propagators is given in Fig. [3l The related FRG equa- 
tions for the propagators are displayed in Fig. [H The 
vertices again satisfy DSE/FRG equations which contain 
vertex functions with a larger number of external legs and 
so on. Thus in practice, solving DSEs/FRGs for general 
momenta requires a truncation scheme. Such truncation 
schemes will be discussed starting with section lTVl Many 
statements in the remainder of this section apply to the 
correct solution, and may be violated by any given trun- 
cation. E. g., any given truncation scheme will in general 
violate constraints, which are imposed by the symme- 
tries of the theory. Furthermore, any truncation scheme 
will neglect part of the physics. The best which can be 
achieved by truncation is that for a given topic these vi- 



olations are sub-leading effects. Depending on the ques- 
tion to be answered and the truncation the irrelevance 
of violations may be established a-priori or a-posteriori. 
This will be discussed further in the following sections. 

In Landau gauge, there is another property of the 
DSE/FRG equations for general amputated Green's 
functions. The equations for fully transversal Green's 



functions v'^Z^^l^.^^^ close on themselves, 



-^(m,n) 
transversal 



FRG/DSE[{F 



transversal 



}]■ 



(10) 



This follows directly from the fact that all internal gluon 
legs are transversal due to the transversality of the gluon 
propagator in Landau gauge, and transversal projections 
on the external gluon legs leaves us with purely transver- 
sal vertices on the rhs of the DSE/FRG equations. It 
is this subset of equations, (dU]), which carries the full 
dynamics of the theory. 

The DSE/FRG equations for Green's functions with at 

least one longitudinal gluon leg, '^'i^aglndmai^ depend both 
on the set of longitudinal as well as transversal Green's 
functions. 



longitudinal 



= FRG/DSE[{F 



{?7l,n) 

longitudinal 



and hence is not closed. The T 



{m,n) 



1 r-pCm,n) |i 
-I ' 1 transversal -J J ' 

(11) 

also have to sat- 



longitudinal 

isfy STIs, that in general imply non-trivial cancellations 
between different diagrams in the DSE/FRG equations 
in pip . This will be be discussed in detail in section [1111 
Here wc discuss standard perturbation theory as a 
prime example for these non-trivial cancellations: In 
leading order perturbation theory, in which only the one- 
loop graphs on the right-hand side in Fig. [3] or the inte- 
grated flow of Fig.dlcontribute, none of the loops is trans- 
verse individually, and only their sum is, up to higher 
order in the coupling constant. 
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D. Analytical solutions and confinement 

There are so far two limits in which an analytical de- 
termination of the leading contribution of Green's func- 
tions can be performed without referring to ad-hoc trun- 
cations. 

One case is the far ultraviolet, where, due to asymp- 
totic freedom, it is possible to determine the leading be- 
havior of Green's functions with analytical means. This 
leading part can be determined with great accuracy, and 
including higher order corrections of the leading part. 
Asymptotic freedom also grants the luxury of proving 
that the leading part of the Green's functions, as de- 
termined by perturbation theory, is in fact uniquely the 
correct one. 

The other kinematical limit is that of the far infrared, 
for external momentum scales 



i.e 



< A^cD- There, a 



general power law behavior of the dressing functions of 
one-particle irreducible Green's functions with 2n exter- 
nal ghost legs and m external gluon legs has been derived 



p(n,m)j-^2-j ^ j-^2-j(n-m)K+(l-n)(d/2-2) 



(12) 



Here, d is the space-time dimension. Below we shall 
restrict ourselves to the most important case d = A. 
Eq. (fT2l) is the unique self-consistent 'scaling' solution 
of the full, untruncatcd tower of DSEs and FRGs [ld| . 

An important consequence of (|12p is the presence of a 
nontrivial infrared fixed point in the running couplings 
related to the primitively divergent vertex functions of 
Yang-Mills theory Q : 



a^^ip') 



An 



gHp')z{p^ 



47r 



COnstgh~gl 

const^g 
const/^g 



(13) 



for p2 ^ 0. Here G{p^) is the dressing function of the 
ghost propagator and Z{p^) the corresponding one for 
the gluon. The appearance of these infrared fixed points 
is independent of the precise value of k, although the nu- 
merical values of the pre-factors are infiuenced indirectly 



with K = K,c = ~kaI'2 in four dimensions. For this solu- 
tion the anomalous dimension k is known to be positive 
0, [3 1 > 0, one therefore finds an infrared divergent 
ghost dressing function and an infrared vanishing gluon 
dressing function. For k = 1/2 the gluon propagator ([14]) 
is finite at zero momentum, < -D(O) < oo, whereas for 
K > 1/2 even the gluon propagator is vanishing in the 
infrared. 

Note that the conservation of global color charge im- 
plies that > in the Kugo-Ojima confinement sce- 
nario [5^. This already enforces scaling [l3| with the 
relation k.a = —'^i^c- In fact the Gribov-Zwanziger sce- 
nario of confinement [ssl . [sgI . Ist} even implies k = kc = 
— ka/2 > 1/2. In addition, a confinement criterion for 
quarks was recently put forward that links the infrared 
behavior of ghost and gluon propagators to the order pa- 
rameter of quark confinement, the Polyakov loop [33|. 
Here quark confinement results in the constraint k > 1/4 
for the scaling solution jS^]. Numerical and analytical 
infrared solutions for propagators obtained in truncation 
schemes [l|, 0, [2^ [H, [13] satisfy the set of constraints 
by the Kugo-Ojima and the Gribov-Zwanziger scenario 
as well as the confinement criterion of [ssj . 

The absence of scaling implies the decoupling of (some) 
degrees of freedom. A solution of this type has been 
discussed in [H, [H, [13, [H, [H, S [13, [H, [sl, 113 and is 
given by: 



Gip" 



const. 



(16) 



We refer to this type of solution as the 'decoupling so- 
lution'. It does not satisfy the set of constraints by the 
Kugo-Ojima and the Gribov-Zwanziger scenario, but sat- 
isfies the quark confinement criterion of [33| : Sk^i— 2k(7 < 
—2. Note that the latter criterion is independent of 
the scenario and comes from the demand of a vanish- 
ing order parameter for quark confinement. We also em- 
phasize that the difference between the scaling solution 
(fT5)) and the decoupling solution is not the possi- 
ble appearance of an infrared finite giuon propagator, 
K.A = —1, but the presence or absence of the scaling re- 
lation KA = — 2kc- We will come back to this point in 
detail in scction lTVl where we discuss numerical solutions 
of the Dyson-Schwinger equations for the ghost and gluon 
propagators. 



In terms of gluon and ghost propagators 
DM = [5,.~^-^^D{p') 

. _ Pj^Py \ Z{p'^) 

n ( ^ ^^^^ 
Dg{p) = — 



(14) 



the solutions yield the power laws 

Z(p2)^(p2)-... G(p2) ^ (p2)-Kc (15) 



III. SLAVNOV- TAYLOR IDENTITIES AND 
RENORMALIZATION 

A. Slavnov- Taylor identities 

We now discuss the relevance of Slavnov- Taylor iden- 
tities when it comes to identifying acceptable solutions 
from functional equations. We present general arguments 
why in transverse gauges these identities cannot serve to 
generate constraints on the transversal Green's functions 
of the theory. This is then demonstrated explicitly for 
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the case of the infrared behavior of the ghost dressing 
function. 

Slavnov- Taylor identities (STI) relate various Green's 
functions. STIs are most conveniently represented in 
terms of the Zinn-Justin equation for the effective action, 
see e.g. fH,!!!, 

f f ST ST ST ST df.Al ST\ 

/ T-—-—. h "T— T 1 T-^T — — cut-off tcrms , 

J ySK^SA'^ SL'^ Sc"' i Sc°- J 

X 

(17) 

together with translation invariance of the anti-ghost, 

X 

Here ^ denotes the parameter of linear covariant gauges, 
which has to be set to zero for Landau gauge. The cut- 
off terms on the right hand side of eq. ([T7]) depend on 
the chosen renormalization scheme and are in particular 
present when a momentum cut-off is used, as necessary 
in all systematic numerical treatments. We come back 
to this point in subsection IIII Bl below. Within the FRG 
approach the rhs of (fT7)l is given in a closed form as a 
one loop term similar to the flow equation see [H, 
Hllljlii,!!^. In general the identities H?]) and ^ are 
solved in a given truncation to the effective action T that 
lead to truncated FRG/DSE-systems for propagators and 
vertices. 

In Landau gauge the constraints from the STIs de- 
couple from the dynamics of the system, if no further 
assumptions are made [g^, HHI- The reason is that the 
STIs relate the longitudinal parts, or more precisely the 
BRST-projccted parts, of vertex fmictions and propaga- 
tors to loops of the transversal gluon propagator, the 
ghost propagator and vertices, schematically written as 

p(m,n) _ orp-rr|-p(m,ri) i r-p()ri,ri) |i .,„n 

longitudinal "-^ It longitudinal/ ' t transversal /J ' K^"^ I 

to be compared with the FRG/DSE-equations for the 
longitudinal Green's functions, (jlip . In turn the 
FRG/DSE-equations for the transversal propagator and 
vertices close on themselves, (fTO)) . as already argued in 
section III Gl Formally, pT|) . can be resolved for 

I 



= 



The last term can be dropped in lowest order perturba- 
tion theory due to the explicit factor g, but this is not 
possible in non-perturbative calculations. 

Nonetheless, the standard truncation of the effective 
action used in the literature within numerical investiga- 
tions of FRG- and DSE systems takes into account full 



rKudinaiPtorsvcrsai]- uou-Abchan gaugc thcorics 
the above hierarchy does not close, there is no self- 
contained subset of relations for a subset of longitudinal 
vertex functions {r|™)J[j„dina,i}j whereas in Abelian theo- 
ries this is indeed possible [63 • 

However, most truncations rely on an expansion in ver- 
tex functions and hence n,m < rimax, JTimax- Then also 
the hierarchy of non- Abelian STIs might be closed. In 
such a case the STIs ([T^ can be resolved and provide 
integral constraints for the transversal Green's functions 

-p()ri,n) 
transversal ' 

We conclude that strictly speaking the STIs can only 
fix the longitudinal parts of IPI correlation functions. 
If, however, the following further regularity assumption 
would be implemented. 

|9pT("'")| < 00, (20) 



transversal and longitudinal correlation functions are 
linked, and the STIs also directly constrain the transver- 
sal correlation functions, for a detailed discussion see [i^. 
Such assumptions have been widely used in Abelian the- 
ories, see e.g. [sl, [H, H^l • In non- Abelian theories results 
obtained from ([20]) have to be discussed cautiously as ([20]) 
certainly fails already at the level of the propagators. 

Having said this, we proceed by exploring the conse- 
quences of P7)) for the correlation functions. Instead 
of using the functional identity ([T7)) that allows for a 
systematic discussion of the truncation scheme involved, 
the STIs for correlation functions are sometimes also de- 
rived by BRST variations of vacuum expectation values 
of product of field operators. 

Similar to functional equations, the STIs relate Green's 
functions with a different number of external legs. We 
would like to discuss this fact within a relevant exam- 
ple. The BRST variation of the expectation value of 
A'^A'^c'^ will produce a Slavnov- Taylor identity involving 
the ghost and gluon propagators, the three-gluon vertex 
and the ghost-gluon vertex, and the irreducible ghost- 
gluon scattering kernel. A simple way to obtain this 
identity is by using BRST invariance on the expectation 
value < A'^A'^c'^ >, yielding 



I 

momentum-dependent propagators, and RG-improved 
ghost-gluon, three-gluon and four-gluon vertices. All 
higher vertices, in particular the ghost-gluon scattering 
kernel, are usually left out. 

As the DSEs for the propagators depend on the bare 
and full ghost-gluon and three- and four gluon vertices. 



- < 4c^9pc" > - < A^^c^d.c^ >+7< ^t^l^p^P > -9ir'^" < A^Alc'^c^ > < A^Alc^c^ >). (21) 
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but do not depend explicitly on the ghost-gluon scatter- 
ing kernel (A^cc), this is apparently a viable truncation. 
Moreover, the latter vertex has been also left out in most 
applications to vertex DSEs, which labels such a trun- 
cation as a minimally self-consistent DSE-approach in a 
vertex truncation. 

However, the scattering kernel does appear in diagrams 
in the FRG equations for both the gluon and ghost propa- 
gators. This comes about as both hierarchies of equations 
represent different resummation schemes and using the 
same vertex truncation at the level of the effective action 



implements different truncation schemes. Note that it is 
precisely this structure that allowed for the uniqueness 
proof of the IR-asymptotics from both sets of equations 
[lO[ , and is also the key ingredient for the important con- 
sistency checks for the results of correlation functions. 

Let us nonetheless proceed and neglect the four-point 
function appearing in (|21J). This has been done previ- 
ously, cither completely [U, |T^, [lO, H^l or only its con- 
nected part |68| . 

Following [67 1, as has been done in general [H. [Tgl. [20| . 
the STI takes the form 



J 



(22) 



r 



Herein F^j/ is implicitly defined by the ghost-gluon vertex gument more. By virtue of Lorentz symmetry it can be 



^u'^\p,q;k) = PuYvt,{p,q;k) 



expanded as [67[ 

Fy^(fc, q\p) = 5^ya - q^^pyh + Pf_,Kc + p„kf_,d + k^ki,e, 



where k is the gluon momentum. It is the ghost-gluon 
scattering kernel, but with two external legs contracted, with five unknown functions a,...,e. 
and thus not the same quantity as appearing in (j21[) . 
which is not contracted and has thus one momentum ar- 



Whcn contracting the left-hand side of with p^kp, 
it vanishes [67!| • From this the relation 



J 



G(g^) 

Z(fc2)' 



■ip, q, k) - 



-a{q,p, k) -pq ^,,^. b{p, q, k) - 6(g,p, k) 



Z(k^) 



\Z{k^) 



Z(k^) 



r 



-pk 



Z(F) 



d-ip, q,k)-qk 



Z{k^) 



d{q,P,k) = 0, 
(23) 



and cyclic permutations thereof, have been obtained [67| . 
However, this condition is actually weaker than possi- 
ble: By the assumption of a tree-level color structure, 
the three-gluon vertex has to be totally antisymmetric 
in any pair of momentum and Lorentz- indices. Hence, it 
is sufficient to contract p2p just with cither pp or fc^ to 
make the left-hand side vanish. Since on the right-hand 
side always one term is transverse w.r.t. the momentum 
which has been used to contract it, it follows that, e. g. 
when contracting with kp that 







Z{p') 



{S^ivP^ - p^iPy)^up{p, q; k)kp. 



This leaves the condition (after relabeling) 

a{p, q, k) — pq b{p, q, k) + pk d{p, q, k) — 0, 



(24) 



and cyclic permutations. This is a condition implying 
certain cancellations between the dressing functions of 
the ghost-gluon vertex, independent of the values of the 
dressing functions of the propagators. In particular, the 



structure is very reminiscent of the one which will be 
obtained in section IIVI for the longitudinal part of the 
ghost-gluon vertex to ensure transversality of the gluon 
propagator in the far infrared. 



Furthermore the equality ([24]) implies the relation (|23p . 
independent of the value of the propagator dressing func- 
tions. The dressing functions a,...,e in eq. (|23p are an (ar- 
bitrary) splitting of the two independent dressing func- 
tions of the ghost-gluon vertex [63], and therefore can 
exhibit in principle any type of behavior in the infrared 
satisfying ([M|) . However, as the non- longitudinal compo- 
nent of the ghost-gluon vertex is infrared finite, while the 
longitudinal one seems to be vanishing Q, the behavior 
of a,...,e is not necessarily regular (as assumed in (20j). 

We conclude that these results do not provide any con- 
straints on the ghost dressing function. Furthermore, the 
condition (j24[) has been derived from the approximation 
([^^ to the complete STI ([^T|) and hence cannot be used 
to exclude certain types of solutions for the Green's func- 
tions as has been attempted in (20j . 
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B. Renormalization 

We now come back to the issue of renornialization- 
dcpcndcnt terms in the STIs. The functional equations 
displayed in Figs. [3] and 2] and the corresponding STIs ap- 
ply to the renormalizcd propagators and hence are equa- 
tions for finite Green's functions. Whereas this finitcness 
trivially carries over to general truncation schemes for the 
functional RG equation due to the regulator insertion, it 
is not evident in general truncations of the DSE. The rea- 
son for the latter intricacy originates in the fact that the 
finiteness of the propagator DSEs is achieved by cancella- 
tions of divergencies stemming from different diagrams in 
Fig. [3l In gauge theories these cancellations arc partially 
governed by Slavnov- Taylor identities. We conclude that 
finiteness of the DSEs requires a careful discussion of the 
renormalization procedure as well as the Slavnov- Taylor 
identities in the presence of a given regularization. 

For any numerical implementation such a regulariza- 
tion involves a momentum cut-off. It is well-known that 
a momentum cut-off scheme might necessitate additive 
renormalization. in particular in gauge theories where 
a momentum cut-off deforms the Slavnov- Taylor iden- 
tities. Indeed, within the flow equation approach the 
momentum cut-off is explicitly introduced at the level 
of the functional integral and the corresponding modi- 
fied Slavnov- Taylor identities are explicitly known [2^ . 
Furthermore it has been shown how to deduce a self- 
consistent renormalization scheme for DSEs with addi- 
tive counter terms in the presence of a momentum cut-off 
from the corresponding FRG equations [13] . In particular 
the latter entail a gluonic mass term that tends towards 
zero if the momentum cut-off is removed. For the DSEs 
this translates into quadratic divergencies that have to 
be canceled identically by mass- type counter-terms (69| . 
This leads to a well-defined unique treatment of trunca- 
tional divergencies in the DSEs that renders the equation 
finite within a given general truncation. In particular 
these considerations enable to interpret the longitudinal 
part of the gluon propagator DSE: It has to be fully can- 
celed by the related longitudinal counter-term in order to 
restore transversality of the fully renormalizcd propaga- 
tor. In turn this uniquely fixes the transversal part. Pre- 
viously, this procedure has implicitly been applied within 
numerical apphcations to DSEs [sl. IbqI. [toI. [80| . 

Having said this we emphasize that multiplicative 
renormalization should still be possible for the DSEs 
within specific schemes. Indeed such a scheme is much 
wanted for as the modified Slavnov- Taylor identities for 
the DSEs in the presence of non-transversal counter 
terms are only implicitly known via the relations to the 
according flow equations. In a multiplicative scheme the 
standard Slavnov- Taylor identities can be used which 
completely decouple the longitudinal parts of the Green's 
functions from the DSEs of the transversal Green func- 
tions and render the latter system fully self-contained. 
Consequently such a choice optimizes the physics con- 
tent of a given truncation. Moreover, for a multiplicative 



scheme it is possible to renormalize the complete system 
with five renormalization constants in Landau gauge, the 
two wave-function renormalization constants of the ghost 
and the gluon, and Z3, as well as the three primi- 
tively divergent vertices, the ghost-gluon vertex with Zi, 
the three-gluon vertex with Z\ , and the four-gluon vertex 
with Z4. The latter three can all be related to a single 
renormalization constant, Zg of the coupling g^. 



IV. GHOSTS AND GLUONS FROM 
DYSON-SCHWINGER EQUATIONS 

In this section we report on explicit solutions for the 
Dyson-Schwinger equations of the ghost and gluon prop- 
agators. These equations involve three fully dressed ver- 
tices; the ghost-gluon, three-gluon and four-gluon ver- 
tices. For the purpose of the present work we shall re- 
strict ourselves to the well-established method of con- 
structing ansaetze for the vertices taking into account as 
much information from the vertex DSEs as possible. We 
would like to emphasize, however, that by now we have 
accumulated enough information on the details of these 
vertex equations [1, 0, [HI, [H, [H, S, [H, [H such that a 
simultaneous treatment of the DSEs for propagators and 
vertices is within reach in the near future. 

Here we investigate the possibility of two different 
types of solutions for the ghost and gluon propagators in 
the infrared. In the course of this study we need to care- 
fully address technical issues like boundary conditions on 
the DSEs, transversality of the gluon propagator in Lan- 
dau gauge and the appearance of quadratic divergencies 
in the gluon DSE once a hard momentum cutoff is in- 
troduced for the numerical treatment of the equations. 
To illustrate our points we use two different truncation 
schemes for the DSEs. The first one has been devel- 
oped in Refs. [l,[23l and produced results in quantitative 
agreement with corresponding lattice calculations in the 
ultraviolet and low-momentum regions. Deviations in the 
mid-momentum regime around one GeV are of the order 
of 20 percent and well understood from the nature of the 
truncation scheme. Qualitative differences in the deep 
infrared are at the heart of our study here and will be 
addressed below. The second scheme will be constructed 
in the next subsections. It serves to study the influence 
of technical issues related to (the breaking of) gauge in- 
variance on the possible types of solutions of the DSEs. 
In turn we will address a connection between renormal- 
ization and boundary conditions in the ghost DSE and 
the issues of quadratic divergencies and transversality in 
the gluon DSE. 



A. A boundary condition in the ghost DSE 

The ghost DSE is given by 
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with the bare ghost-gluon vertex rj^^'"'^''^''^ its dressed 
counterpart F/j ' and the momentum routing k — p + q. 
The ghost renormaUzation function and the ghost- 
gluon vertex renormaUzation function Zi both depend 
on the renormaUzation point 

The equation as it stands is multiphcatively renormal- 
izable as recaUed cxpUcitly in appendix|B] It then foUows 
that the choice of g^{iJ?), which impUcitly fixes can- 
not have any impact on the type of solutions one finds in 
the DSEs. All multiphcatively renormalizable solutions 
must exist between infinitesimally small (asymptotic 
freedom; large /i^) and a potential maximum value that 
can be derived from the infrared behavior of the cou- 
plings (jl3p . To discuss this further wc rewrite the ghost 
DSE schematically as 

G(^)=^'-^'^'^^'^^^^P'^^'^ (26) 
and subtract the equation at = 

G^) - G(0i7)-^^^'^'^'^^'^^='^^~'^'^'^^^- ^''^ 

This way we got rid of the renormaUzation factor Z3 in 
exchange for the boundary condition 0(0; /i^). One can 
now explicitly choose between solutions of the scaling 
ty pe (1151) or the decoupling type pB]) by varying G(0; /i^) 
[J, [2^. Clearly, G(0;/^^)~^ — corresponds to an in- 
frared diverging ghost dressing function. This choice cor- 
responds to the 'horizon condition' derived by Zwanziger 
to account for the presence of the Gribov-horizon when 
evaluating correlation functions [H, [H, [7l| . The other 
option, G(0;//^) = const., produces an infrared finite 
ghost dressing function by construction. If the un- 
subtracted equation (j26p is used these conditions trans- 
late into equivalent ones for the renormaUzation constant 
Z^. From this it is clear that it is instructive but not 
necessary to subtract the ghost equation exactly at zero 
momentum. One can find both solutions for an arbi- 
trary subtraction point s, if the boundary conditions on 
G(s; /i^) arc selected appropriately. This is a consequence 
of renormaUzation group invariance. For the scaling solu- 
tion the equivalence of boundary conditions for G{0\jj?) 
and G{s;fj,^) has been explicitly demonstrated in |72l |. 

In [2^ the two types of solutions have been shown 
to coexist in a numerical treatment of the ghost DSE 
alone. Below we demonstrate that self-consistent solu- 
tions of both types exist when the ghost and gluon Dyson- 
Schwingcr equations are solved simultaneously. 



I 

B. The gluon DSE: quadratic divergencies 

Having discussed the ghost DSE on general grounds 
we now focus on the gluon DSE. To this end we need to 
specify an explicit truncation scheme for the fully dressed 
vertices involved. As stated above, one such scheme has 
been constructed in refs. [i,[23. We first shortly summa- 
rize this truncation before we qualitatively improve upon 
it. 

(2 1) 

In general, the dressed ghost-gluon vertex FJi ' can be 
written as 

F(,2.1) (p, q; k) = iq,A{q. k) + ik^B{q, k) (28) 

where p and q are the incoming and outgoing ghost mo- 
menta respectively and k is the gluon momentum. A{q, k) 
and B{q, k) denote the two dressing functions of the ver- 
tex with A{q, fc) — *■ 1 and B{q, fc) return the tree- 
level vertex. The dressing functions A(g, k) and i3(g, k) 
have been subject to investigations in the continuum 
[1,0] and A{q, k) also on the lattice [5l|, [s^l and found to 
only mildly deviate from the tree-level behaviors. Fur- 
thermore it has been shown in Q that the infrared be- 
havior of the ghost-gluon system is only mildly affected 
by possible vertex dressings. This justifies the truncation 

Tf^\p,q;k)^iq^, (29) 

i. e., the replacement of the fully dressed ghost-gluon 
vertex in the DSEs by its tree-level counterpart. This 
truncation has been used in 0, [t^I and [l^l • 

The situation is somewhat more complicated for the 
three-gluon vertex. Here an ansatz v'f^^'^p {j),q) = 
F^°i^p''^°''r^3(p, q) has been chosen in Q with the tree- level 
vertex F^'l/p'*'''-' and dressing function F'^^(p, q) that leads 
to solutions for the ghost and gluon DSE with the correct 
one-loop running from rcsummcd perturbation theory. 
The ansatz is given by eq. IAl|in the appendix. Further- 
more it has been shown in [70[ that this choice of the 
vertex also leads to solutions which respect multiplica- 
tive renormalizability of the gluon DSE. Contributions 
involving the four-gluon interaction have been neglected 
in [H. For the scaling solution this affects only the in- 
termediate momentum regime, as shown there. For the 
decoupling solution the omission of the four-gluon vertex 
may also have an impact at low momenta. This needs to 
be explored in future work. 

In general a truncation scheme such as the one sum- 
marized here faces the problem of quadratic divergencies 
appearing in the gluon DSE. These have been dealt with 
in the past by explicit subtraction procedures on the level 
of the integrands 0, [tI, HO] , numerical subtractions on 
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the level of the right hand side of the gluon DSE [69| or 
subtractions of purely quadratically divergent term in- 
volving only the perturbative form of the ghost and gluon 
dressing functions [T3|. Although conceptually different 
all these procedures lead to similar results if correctly 
implemented. 

Here we introduce yet another procedure which takes 
into account the general considerations of subsection 
IIIIBI Recall that the introduction of a hard momen- 
tum cut-off into the theory modifies the Slavnov- Taylor 
identities. This modification has its biggest effects on 



momentum scales near the cut-off. An improved trunca- 
tion scheme that avoids the related problem of quadratic 
divergencies therefore should work with vertex modifica- 
tions at these scales. Such a scheme is constructed in the 
following. We devise ansaetze for the ghost-gluon and the 
three-gluon vertices that resolve the issue of quadratic di- 
vergencies. The resulting ghost and gluon DSEs arc then 
multiplicativcly renormalizablc in the first place without 
the need for additional subtraction procedures. 

The procedure turns out to be straightforward. For 
the ghost-gluon vertex we start with the general ansatz 



/ 2 2 2 

Tf^\p,q,k) ^ ^qJl + A + B^ + C^ + D^ + 



which includes all possible modifications of the vertex in 
terms of dimensionless combinations up to powers of two 
in momenta and angles. Here p is the incoming ghost mo- 
mentum, q is the one of the outgoing ghost and k ^ p — q 
is the momentum attached to the gluon leg. The coef- 
ficients A.. J are assumed to be constant in momentum. 
One can then perform an ultraviolet analysis of t he g luon 
DSE and the ghost DSE along the hues of ref. [1,[70| and 
restrict the possible values of the coefficients by demand- 
ing the following conditions: 

• No quartic and quadratic divergencies should ap- 
pear in the ghost loop of the gluon DSE and the 
ghost-gluon loop in the ghost DSE. 

• The logarithmic divergencies of these loops should 



pz pz qz p'^q-^ p^k'^ q'^k'^ J 

(30) 

I 

remain the same as for the bare ghost-gluon vertex, 
i.e. as in perturbation theory. 

This restricts the coefficients of (pO|) to a family of solu- 
tions which depend on the explicit values of two of the 
constants A.., J. Setting as many of these constants to 
zero as possible one arrives at the surprisingly simple 
form 

f(y)'^^b,g,fc) = zg,(l-^) . (31) 

A similar procedure for the thrce-gluon vertex produces 
the result 



f(°i)'^^(,,p)==r(L^)(°)(5,p)r^^(p,, 



140 



52 p2 
T7P 



sy p^ 
51 g2 



52 q2 
17fc2 



26 P 
17 g2 



104 {q.kf 

17 92fc2 



(32) 



where Tp^J {q,p) is the bare three-gluon vertex and 
T'^^ jp. q), given in appendix \X[ has been introduced in 
[to] to reproduce the correct logarithmic running of 
the gluon in the ultraviolet. The additional, power- like 
vertex dressings in pip and (|32p serve to avoid quadratic 
divergencies in the ultraviolet momentum regime. As 
discussed above, the power-like modifications should af- 
fect the vertices only for momenta close to the ultra- 
violet regulator. We therefore multiply ([3T|) and ([32|) 
with appropriate damping factors that cancel these terms 
smoothly with decreasing momenta. The explicit form of 
the damping factors will be given in the next subsection. 



C. Transversality 



Next we come back to the issue of transversality, con- 
tinuing the discussion from section IIII Al In Landau 
gauge, the gluon propagator is transverse and this prop- 
erty is reflected by the structure of the full, untrun- 
cated gluon DSE. For symmetry conserving regulariza- 
tion schemes this property can be maintained in trun- 
cation schemes that satisfy the Slavnov- Taylor identi- 
ties. However, even truncation schemes satisfying a finite 
number of STIs are extremely hard to implement in nu- 
merical calculations. There one usually has to resort to 



11 



a momentum cut-ofF. As a consequence spurious longi- 
tudinal terms are generated in the gluon DSE that need 
to be projected out by contracting the gluon DSE with a 
transverse projector. This problem has been overlooked 
in [l3|. There a truncation has been used that is for- 
mally transverse under dimensional regularisation while 
it is not when a hard cutoff is used in the numerical cal- 
culations. Consequently, in contradiction to their claims, 
the numerical solutions presented in (l7j are not strictly 
transversal. 

In general, projection onto the transverse part of the 
right hand side of the gluon DSE, as done in 0] guar- 
antees reliable results. Nevertheless it is interesting to 
construct a truncation scheme that in addition minimizes 
the spurious longitudinal terms. This is the purpose of 
this subsection. 

As is apparent from the analysis of ref. 0] transver- 
sality is naturally maintained in the ultraviolet where 
perturbation theory is at work. In the infrared, how- 
ever, the situation is different and sizeable longitudinal 
components of the gluon are present. To improve this sit- 
uation we consider the most general form for the dressed 
ghost-gluon vertex, eq. , and contract the ghost-loop 
in the gluon DSE with a longitudinal projector, thereby 
projecting onto the spurious terms. Apart from trivial 
factors one obtains the combination 

q,^T(^''Kp,q;k) = q/-lt^{k,A{k,p)+p,B{Kp)) 
p'^ p^ 

== (33) 

where p is the gluon momentum and q and k = q ~ p 
are the two ghost momenta. This is reminiscent of 



equation p4|) . From this relation we find the condi- 
tion B{k,p) = —A(k,p)^ which eliminates the spurious 
terms and makes the ghost-loop transverse. As will be 
shown later this construction is sufficient to render the 
gluon DSE transversal on the level of numerical accuracy. 
Note, however, that a more general construction for B is 
possible such that exact transversality of the gluon DSE 
for all momenta is guaranteed. 

As a result the dressed ghost-gluon vertex in the in- 
frared can be written as 

f 1,2,1) q. k) = ^q^A{q, k) - ikJ^A{q, k) (34) 

in the momentum labeling of eq. (j28p where k is the gluon 
momentum. It is important to note that the additional 
term in eq. p4p is longitudinal in the gluon momentum, 
i.e. it does not contribute to the ghost-loop in Landau 
gauge but merely serves to eliminate non-transverse arte- 
facts. Thus the infrared analysis of the ghost loop per- 
formed in [J] and the value of k remain unchanged. 

Longitudinal contributions to the ghost-gluon vertex 
have been calculated from the vertex-DSE in 0] and 
found to vanish for large momenta, in agreement with 
perturbation theory. We therefore need to multiply the 
longitudinal structure in (|34p by appropriate damping 
factors to account for this behavior. Furthermore we 
choose A{q, fc) = 1 as in the truncation scheme of the 
previous section. 

Taking into account the ultraviolet modifications (|3ip 
and (j32p of the previous subsection we then find our final 
ansaetze for the ghost-gluon and three-gluon vertices 



m\q,v) 



fuvip,q;k) -ik^ 



r(0'3)(o)(^,p)r3.(,,p) 



k.q 



fc2 

140 52 p2 



fiR{p,q;k) 



89 p2 
51 g2 



52 g2 



26 P 
17^ 



104 {q.kf 
17 q^k^ 



(35) 



fuv{p,q;k) (36) 



with appropriate damping factors 



fuvip,q;k) = 



(p' 



A^ 



)iq' 



fiR{p,q;k) 



(p2 



A2 

^^IR 



){q' 



(37) 



for the infrared and ultraviolet modifications of the ver- resulting forms of the equations for the ghost and gluon 
tices. The dressing F'^f (p, q) is given in appendix [51 The DSEs are given by 
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Z(p2) ^ ' " 3 
"3 



Zs + g ^ TTTTI G{q )G{{p-q) ) 



1 - 2, 



G(p2) ^ ^ '7 (27r)4 p'^q 



(27r)4 p2q2 




^ A Qip^qMp- 




(27r)4 p2^2 









/ "^"'^r; ' Z{q')Z{{p~qf)T'<^{q,p) , (38) 



^3-5^iV, ^"'i;:, ^ G(y)Z((p-g)^) . (39) 



These forms are similar to the ones of the previously 
used truncation scheme of fs, 'to'I . The effect of the ver- 
tex modifications ([55]) and ([55]) is hidden in the kernels 
K, M, Q which are given in appendix [X] together with 
the original kernels K, M, Q of ref . Q • 

The two free scale parameters Am and Auv are fixed 
by the following procedure: Once a numerical solution of 
the DSEs has been obtained we project the right hand 
side of the gluon DSE onto the longitudinal component 
of the gluon propagator. This component is then min- 
imized by changing the two scales. We find that this 
procedure is insensitive to Auy in a window of roughly 
A/2 < A(jY < A, where A is the ultraviolet cutoff in our 
numerical procedure. In practice we choose the lower 
bound of this interval. The infrared modification of the 
ghost-gluon vertex is important in the regime where the 
infrared asymptotics sets in. Consequently we find an 
optimized value of Am « 300 MeV. This choice leads 
to spurious longitudinal components of the gluon propa- 
gator in the mid-momentum regime of the order of less 
than one permille, i.e. of the order of the numerical error 
of our calculations. In the infrared and ultraviolet mo- 
mentum regime these terms are even smaller by orders 
of magnitude. 

Finally we wish to point out that the technical im- 
provements of this section are also extremely useful when 
it comes to calculations at finite temperature, where 
quadratic divergencies represent a much more severe 
problem than at zero temperatures. Improvements along 
the lines suggested here have been anticipated in [t^ . 



D. Renormalization conditions 

Having specified the truncation scheme it remains to 
detail the renormalization conditions imposed on our sys- 
tem of equations. To discuss these we need to recall the 
definition of the running coupling from the ghost-gluon 
vertex [l[ 



a{p ) ^a{fj, ) 



2,G(p2;/i2)2^(p2.^2) 



G(Ai2;/i2)2^(^2.^2) 



(40) 



where the dependence of the dressing functions on the 
renormalization point is given explicitly. Note that 
the right hand side of the equations is independent of fi^, 
i.e.; it is a renormalization group invariant. This defi- 



nition follows straightforwardly from the Slavnov- Taylor 
identity 



Zi 



ZgZ^Z^ 



(41) 



where the renormalization function Zi of the ghost-gluon 
vertex is always finite in Landau gauge. For a bare ghost- 
gluon vertex, as chosen here, Zi ^ 1. Then, in the mo- 
mentum subtraction scheme {MOM) defined in [l| the 
renormalization condition G{fP; fP)^Z(fi^; fi^) = 1 is im- 
posed, which we also adopt here. Then G{ij?';ijP') and 
Z{ij?;ii^) are not independent of each other. In practice 
one chooses a value for Q:(/i^), which in turn determines 
the renormalization point . In [23| the special choice 
G(/x2; /^2) — Z{p?; ^2) = 1 has been used, which in gen- 
eral leads to Zi = const 7^ 1 by virtue of the STI (|4ip . 

The definition of the running coupling represents 
a mass independent renormalization scheme. This, how- 
ever, may not be adequate for the decoupling solution 
with < D{0) < 00. In this case one should separate 
the 'massive' part of the gluon propagator from the part 
proportional to p^ . Therefore we introduce 



D-\p') 



P 



1 



Zip') Z{p^) 



(42) 



with the renormalization group invariant mass parameter 
TO. Note that this splitting is not unique and one obtains 
an additional free normalization parameter Z[Q). The 
running of Z then enters the running coupling by 



a[p ) = a{fi )- 



'G(a.2. ^2)2^(^2. ^2 



(43) 



and the renormalization condition is given by 
G{^j?] ^'')'^Z{^'^; ^'^) ~ 1. Independently of this 
renormalization condition one has the freedom to vary 
Z(0), which together with the boundary condition 
G(0) determines the value of the infrared fixed point 
of this coupling. This is further discussed in the next 
subsection. 

For all our numerical results except the ones presented 
in the appendix [B] where we show the independence on 
this choice, we have used q;(/^2) = 1, 
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p" [GeV^] 

FIG. 5: Numerical solutions for the ghost and gluon dress- 
ing function with different boundary conditions G(0). The 
(artificial) longitudinal components of the gluon propagator 
are not displayed, since they are of the order of less than one 
permille, i.e. of the order of the numerical error of our cal- 
culations. All results shown here are obtained from our novel 
truncation scheme. Differences to the scheme defined in [sI-ItoI] 
are, however, only very small and would not be visible in the 
plots. 



E. Numerical results 

Our numerical solutions for the ghost and gluon dress- 
ing functions arc shown in Fig. [S] The corresponding 
momentum scale has been fixed by best-possible match- 
ing of the gluon dressing function to the corresponding 
one on the lattice, cf. 



V. SEC:LATTICE 

. Thus we inherit the lattice scale. All results displayed 
are obtained from our novel truncation scheme. Differ- 
ences to the scheme defined in [t^I are, however, only 
very small and would not be visible in the plots. This pro- 
vides additional justification that the old scheme already 
represented a reliable result. Trans versality is manifest in 
our new truncation scheme; the longitudinal components 




[GeV^l 



FIG. 6: The running coupling for the two types of solu- 
tions, defined as a(p^) = a{^^)G(p'^)^ Z{p^) (top diagram) 
and a(p^) — a{fi^)G{p^)^ Z (p^) (bottom diagram). 



of the propagator (not shown in the figure) are smaller 
than one permille and therefore of the same size as the 
numerical error of our calculation. Quadratic divergen- 
cies do not appear and the numerical solutions respect 
multiplicative renormalizability, see appendix [BJ 

By varying the boundary condition G(0,/i^) we are 
able to generate both, a solution of the scaling type I, eq. 
din])), with G"^(0, //^) = and a continuous set of decou- 
pling solutions characterized by a finite value 0(0, /i^). 
Shown are results for three different values of G(0,/i^). 
The corresponding gluon propagator is either massive in 
the sense that D{Q) — \imp2 Z (p^) / = const, for 
decoupling, or has the power like behavior psp with 
K = Kc = {93- yT2()T)/98 « 0.595353 H in the case 
of scaling. In the ultraviolet momentum region, both 
types of solutions are almost identical, as expected. The 
running couplings are shown in Fig. [S] For the scal- 
ing solution one observes the infrared fixed point known 
from previous studies, whereas the coupling for the de- 
coupling solution falls with p'^ in the infrared when the 
mass-independent definition (j40p is used (top diagram 
of Fig. ini). If, however, we employ the mass-dependent 
definition (bottom diagram of Fig. (5]) the resulting 
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coupling develops a fixed point in the infrared. The value 
of this fixed point is an additional free parameter and has 
been fixed at a(0) w 3, i.e. at approximately the value 
of the fixed point of the scaling solution. Note that these 
couplings still show variations with changes of G(0) at 
intermediate momenta from 0.1 — 1 GeV. 

In the light of these findings we wish to make the fol- 
lowing comments: 

• In principle one has to face the problem that an in- 
adequate ultraviolet rcnormalization may result in 
an infrared mass for the gluon reflecting the break- 
ing of BRST-invariance by the momentum cutoff. 
We explicity checked that this is not the case for 
both our scaling and the decoupling solution; the 
solutions are independent of the cutoff A imposed 
in the loop integrals. In particular this is true for 
D{0). 

• We also assessed the behavior of the dressing 
functions under a change of the renormalization 
point i'^, i.e. we chose Oi{iy'^) = 0.5 
and 0(1/^; i/-^) = 1. Again, by varying 
G(0, v^) we find both types of solutions and the 
resulting ghost and gluon dressing functions are re- 
lated to the ones at a{n'^) = 1 by renormalization 
factors, as they should. The explicit results can 
be found in appendix |B] Wc therefore conclude 
that there is no notion of a 'critical' coupling dis- 
tinguishing the two types of solutions as claimed 
in [201 . Instead it is solely the boimdary condition 
G{0,fi^) that matters. 

• We wish to point out that the decoupling solution 
with an infrared finite gluon propagator necessar- 
ily entails an infrared finite ghost dressing function 
(up to logarithms) as long as the ghost-gluon vertex 
does not contain any nontrivial power-laws [lol . [20| . 
This situation is also realized in [l3| , where hence- 
forth Kc = must be realized. This is compatible 
with part of the fits they deduce from their numer- 
ical solutions. 

• The issue of transversality is separate from the 
question of which types of solutions are realized in 
a particular truncation scheme. In non-transverse 
truncation schemes longitudinal components of the 
gluon propagator only constitute a problem if they 
are back-fed into the equation by the use of a non- 
transverse projection method of the gluon DSE. 
This has carefully been avoided in the truncation 
scheme of Q which consequently delivered very 
similar results as the new scheme employed in this 
section. 

• We wish to emphasize again that either (i) global 
BRST symmetry is unbroken, then the decoupling 
solution would imply the breaking of global color 
symmetry as indicative for a Higgs phase of the the- 
ory. Or, (ii), in a confining phase the decoupling 



solution implies the breaking of global BRST sym- 
metry and therefore does not agree with the Kugo- 
Ojima confinement scenario. Indeed, all known 
BRST-quantizations indicating towards an infrared 
finite ghost even break off-shell BRST [H, SO]. 
Therefore it is not clear how to construct a physi- 
cal state space in the decoupling case. In terms of 
Green's functions this implies that it is not known 
how to construct physical observables in general. 
In our opinion this is a serious problem. In addi- 
tion the decoupling solution is in contradiction with 
the Gribov-Zwanzigcr confinement scenario due to 
its finite ghost dressing function [t^. All in all, 
the status of the decoupling solutions is therefore 
clearly different from the scaling solution, which 
agrees with both, the Kugo-Ojima and the Gribov- 
Zwanziger scenarios. 

It is furthermore interesting to compare our results to 
corresponding ones from DSEs on a torus. In 0,[7H| solu- 
tions from the ghost and gluon DSE on a torus have been 
discussed which did not connect to a scaling type of solu- 
tion even if very large volumes were taken. In a later work 
[ttI ] the renormalization conditions for the torus solutions 
have been reconsidered and adapted such that the scal- 
ing solution in the infinite volume/continuum limit has 
been reproduced. The status of the solutions found in 
0, [zB], however, remained somewhat unclear. From the 
results of this work we arc now in a position to clarify 
this issue. The different renormalization conditions on 
the torus employed in 0, [t^ and 77 1 are in one-to-one 
correspondence to the boundary condition G~^(0, ij?) for 
the continuum DSEs investigated in this work. The in- 
finite volume/continuum limit of the solutions found in 
0, [7^ is therefore given by one of the decoupling type of 
solutions reported in this work. 

Finally we wish to point out that the parameter 
G^^(0, corresponds to a related one in the Zwanziger- 
Lagrangian approach (t^I to infrared Yang-Mills theory. 
There, a mass parameter can be introduced which pro- 
duces decoupling solutions for non- vanishing mass 40 1 
and scaling solutions for vanishing mass, see e.g. |79 | 
and refs. therein. The latter is in agreement with the 
original Gribov-Zwanziger scenario. 



VI. GHOSTS AND GLUONS FROM THE 
FUNCTIONAL RENORMALIZATION GROUP 

The DSE-analysis of the previous sections can be re- 
peated within the functional RG. In order to facilitate 
the access to the FRG-literature on Landau gauge Yang- 
Mills, e.g. [13, m, M, M, H M, mMkM, we 

present the results in standard FRG-notation and pro- 
vide the dictionary to the DSE-computations. The re- 
sults shown here are preliminary results taken from [6^ . 
We shall be brief on the details of this analysis and only 
discuss the relation between the renormalization proce- 
dure in the DSE of section ITV] and that in the FRG. 
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Furthermore we comment on the truncation and discuss 
its relation to that used in the DSE. The renormahzation 
conditions are contained in appropriately chosen initial 
conditions for the effective action Fa at the ultraviolet 
scale A. Indeed it has been shown in [U that the inte- 
grated flow equations define a set of DSEs within a con- 
sistent BPHZ-type renormahzation scheme. This system 
reads for the effective action 



F.FA + lTrlnF(^)-l/fTrFi^ 



.(2) 



Rk' 



(44) 



The first term on the rhs of ([44]) comprises the renor- 
mahzation conditions as an initial condition as can be 
seen via the comparison of the second field derivatives of 
dH]) and (123). The second term on the rhs of (gH) con- 
tains one loop contributions to Green's functions p*^^"'™' 
with full propagators and vertices, the third integral term 
contains RG- improvement terms such as the two-loop di- 
agrams in the DSEs. In the present truncation we need 



the (inverse) propagators F)^ = F*^"--^'. F 



.(2) 



2\„2 cab 



ab 



2 \'2 cab 



= F(2.o), 
longitud. , 
(45) 



as well as the vertices F^^'i) , T^°-^^) , F^"-^) . For the vertices 
we refer to the parameterization in section IIVI and to 
[13, [Hj [131 ■ The wave function renormahzation functions 
ZatZc relate to the dressing functions Z{p'^),G{p'^) as 



Za{p') 



1 



Z(p2) ' 



ZciP^) 



1 



G(p2) 



(46) 



(2.0) 



(0.2) 



and 



whereas Z3, Z3 and Zi are defined with F^^ 

(2 1) 

F\ . Note that RG-invariancc of the FRG-equation is 
automatically guaranteed by the separate RG-invariance 
of each term in (|44)) . and hence the invariancc of the 
FRG solution under a change of the renormahzation scale 
^. The scaling and decouphng solutions (at k = 0) 
are now adjusted via fine-tuning of Z3 as in the DSEs. 
Indeed, the ghost FRG can be mapped into the corre- 
sponding DSE and one can directly use eq. (P7)) . Fi- 
nally, the truncation is optimized 0, by an appro- 
priate choice of the regulator function R. The optimized 
regulator employed for the present work has been de- 
rived from functional optimization in It has been 
already shown for the present truncation in [2^ that 
optimized regulators lead to the same analytic expres- 
sions for k's and a(0) in the scaling solution as for the 
DSE equations, kc ~ 0.595353, a(0) « 3. Other regu- 
lators lead to scaling solutions with slightly varying k, 
K e [.539... , .59535...], for a specific choice see also j33 |. 
Note that the above interval serves as an estimate on the 
systematic error of the present FRG and also DSE trun- 
cation. Furthermore it has been argued in (30j that the 
lower value relates to finite volume studies, providing a 



link to the torus DSE analyses in [i,[2i,[23l- More details 
on the present FRG analysis will be presented in [6^ . 




FIG. 7: Numerical solutions for the ghost and gluon dressing 
function with two different boundary conditions G(0) calcu- 
lated from the FRGs. 
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FIG. 8: The running coupling from the FRGs, defined as 

The numerical solutions from the FRGs are presented 
in Figs. [7] and (5] We show the scaling solution and 



16 



one representative for a decoupling solution. They agree 
qualitatively and to a large extent also quantitatively 
with the ones from the DSEs presented in the previ- 
ous section. Therefore both frameworks support the con- 
clusions discussed there. Quantitative differences in the 
mid-momentum regime can be attributed to the differ- 
ent diagrammatics used in both approaches. In summary 
we have discussed the equivalence and consistency of the 
renormalization procedure for both, DSEs and FRGs. 
Moreover, the FRG provides a consistent momentum cut- 
off regularization of the corresponding DSE equation via 
and thus allows to deduce the modified STIs for the 
DSE in the presence of an ultraviolet momentum cut-off, 
see [13, US] . A crucial difference in the present truncation 
is the tadpole diagram in the gluon FRG-equation that 
depends on the full four-gluon vertex. This incorporates 
two-loop contributions of the sunset diagram in the giuon 
DSE, see Fig.H 



VII. COMPARISON WITH LATTICE RESULTS 

In the previous two sections we obtained two different 
types of solutions for the ghost and gluon propagators in 
the DSE and FRG approaches. It is certainly instructive 
to compare these results to the ones from lattice calcu- 
lations. As became apparent from a number of works in 
the past years such a comparison is not unambiguous. 
Ideally one strives for a situation where exactly the same 
quantities are calculated in the continuum and on the lat- 
tice. However, this is currently not the case for a number 
of reasons. First, lattice calculations are necessarily done 
in a finite volume. It is therefore mandatory to take into 
account finite volume effects and zero mode contributions 
absent in the infinite volume/continuum limit. Second, 
one encounters finite size contributions due to the non- 
vanishing lattice spacing. Third, artefacts due to the 
gauge fixing procedure are different from the ones in a 
continuum formulation. 

Before we discuss these issues further let us com- 
pare the continuum solutions with the lattice results of 
refs. 0, [sij in minimal Landau gauge. In the top dia- 
gram of fig.[9]we display the gluon dressing function from 
different approaches. At large momenta, where pertur- 
bation theory sets in, all results are in excellent agree- 
ment with each other. The DSE results as well as the 
FRG results in the intermediate regime show only a mild 
dependence of the type of solution, i.e. scaHng or de- 
coupling docs not really matter here, as expected. As 
compared to the standard DSE results the dressing func- 
tion from the functional RG approach is closer to the 
lattice data. From the discussion of the last section this 
was to be expected, since the FRG truncation included 
effects from the gluonic two-loop diagrams neglected in 
the DSE-truncation. Note that such contributions can be 
either included directly or phenomenologically by modi- 
fying the three-gluon interaction in the one-loop diagram 
also into the DSE framework, see e.g. (s^ . 
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FIG. 9: Both type of solutions of sections llVl and lVIIc ompared 
to lattice results in minimal Landau gauge from |4J, [81(. 



The infrared behavior of the propagator functions for 
the gluon, D{p^) = Z{p'^)/p'^ , of both solutions are com- 
pared in the second panel of fig. O Clearly, the scal- 
ing solution comprises an infrared vanishing propaga- 
tor, whereas the decoupling solutions are infrared finite. 
Changing the boundary condition G~^(0, /x^) from zero 
to finite values first leads to a finite but small value for 
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D{Q) with the corresponding gluon propagator still be- 
ing non-monotonous. From a certain minimal value of 
6*^^(0, /i^) on, this behavior changes and the gluon be- 
comes a monotonously decreasing function of momen- 
tum. Such a monotonous behavior is also seen in the 
lattice data, which therefore clearly represent a decou- 
pling type of solution for the gluon. 

The third diagram of fig. [5] compares results for the 
ghost dressing function from the continuum approaches 
with the lattice data. Again we find that the lattice re- 
sults resemble a decoupling type of solution, i.e. it seems 
from the plot that the corresponding value 0(0, /i^) on 
the lattice would be finite. 

In conclusion, the lattice results support a decoupling 
solution with broken BRST symmetry. The origin of this 
BRST breaking is currently not understood. It may be 
induced by one or more of the three aforementioned ef- 
fects. 

Finite volume corrections are unambiguously present 
and are sizeable even for large volumes |77i . However, 
recent calculations on very large lattices [46l. l47l. [49j have 
demonstrated that these corrections alone do not change 
the solution from a decoupling one to a scaling one. 
These results have been obtained on rather coarse lat- 
tices in the minimal Landau gauge^. Correcting for both 
effects has a significant impact. 

Recently, it has been shown [sJl that the gluon prop- 
agator, and in particular its value at zero momentum 
Z?(0), are affected by discretization artifacts which orig- 
inate from different possibilities to define the gauge field 
on the lattice. The consequences of such discretization 
artifacts are not yet studied in detail, but evidence has 
been found that is in favor of a scaling solution on very 
fine, large lattices (53 |. 

Finally, the question of adequate gauge-fixing is of rel- 
evance. As discussed in the beginning the basic scenario 
is established in the absolute Landau gauge, and there- 
fore the lattice calculations should be done in this gauge. 
However, the numerical implementation of this gauge is 
expensive, and hence most calculations so far have been 
done in the minimal Landau gauge. Nonetheless, it is 
known that at least the ghost propagator differs in both 
gauges, but the extent and volume-dependence is not 
known sufficiently well in four dimensions (48i . [sij . Again, 
recently [11] evidence has been found which suggest that 
the absolute Landau gauge produces the scaling solution. 
However, only in three and two dimensions support for 
this has been found yet and only indications are so far 
available in four dimensions [H, [H, [s^ , mainly due to 
the computing time required in four dimensions. Extrap- 
olations in the number of dimensions d make it nonethe- 
less likely that also in four dimensions in the absolute 



Note that minimal Landau gauge may have ambiguities, that 
possibly lead to differing propagators or at least differing finite 
volume effects [56l |. in contrast to the unique definition of the 
absolute Landau gauge. 



Landau gauge the propagators are scaling. 

Hence both effects are in favor of a scaling solution. 
Nonetheless, the interplay of these effects and finite vol- 
ume effects has not yet been fully investigated, and the 
situation in lattice gauge theory is not yet unambiguously 
resolved. 

By comparison of the lattice and the functional re- 
sults it seems furthermore tempting to conjecture that 
the scaling solution is in fact obtained in absolute Landau 
gauge, and the absolute Landau gauge is implemented 
in functional calculations by the implementation of the 
boundary condition G(0, ~ 0. The good agreement 

of the massive solution to the lattice results in minimal 
Landau gauge, see Fig. [9l suggest further that the other 
boundary condition G{0, y.^)^^ 7^ < oo seems to imple- 
ment the minimal Landau gauge. However, final justifi- 
cation of this idea requires a clarification of the role of 
discretization artifacts [13] and other subtleties in defin- 
ing the minimal Landau gauge [56| . 

Alternatively, one could check which solution satisfies 
the condition for absolute Landau gauge. One is the pos- 
itivity of the ghost dressing function, which is given in 
both cases. The other is the requirement of minimiz- 
ing the trace of the gluon propagator, equivalent to the 
condition dH) [H], 

T{D) = c J pZ{p)dp, (47) 

with a positive constant c [11] . As is visible from Fig. \E\ 
and [21 at low momenta the scaling solution is smaller 
for both the DSE and FRG results. At mid-momentum, 
both solutions are similar from the FRG, while from the 
DSEs the massive solution is somewhat smaller. This 
already shows that the question whether (|47p is minimal 
for scaling or decoupling depends on the details of the 
truncation and therefore cannot be answered from the 
results plotted in Fig. [5] and [HI 

Finally, we would like to mention two general issues. 
First, we expect that the general structure of the inter- 
play between gauge fixing and the infrared behavior of 
Green's functions is general and also present in other 
gauges. Clear indications for such a picture are given 
in [85| for the case of Coulomb gauge. Second, we wish 
to emphasize that our discussion is independent of the 
number of colors, Nc, and therefore valid for SU(7Vc)- 
Yang-Mills theory. In the DSE/FRG framework an in- 
direct and very mild dependence of the ghost and gluon 
dressing functions on Nc is induce d by sub-leading com- 
ponents of the four-gluon vertex |12| . which only have 
an impact in the mid-momentum region. This mild de- 
pendence is also found in lattice calculations comparing 
the dressing functions for the cases of SU(2) and SU(3) 

MM- 

In fact, based on the structure of the DSEs it has 
been conjectured [11] that the scenario presented here 
can be extended to any (scmi-)simple Lie-group without 
any qualitative change. First support for this conjecture 
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in the case of the exceptional Lie-group G2 in lattice cal- 
culations has been presented in [89|. 



VIII. GLUON CONFINEMENT 
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FIG. 10: The absolute value of the Schwinger function A{t) 
plotted against time for both, the decoupling and scaling type 
of solutions. (The latter result with a slightly different scale 
has previously been published in [90l].l 

Finally we wish to investigate the issue of positivity 
violations in our two types of solutions. The confinement 
of gluons by the absence of a spectral representation of 
its propagator has been addressed e.g. in [l|, [7l|, [9^ with 
violation of positivity being a sufficient criterion for the 
absence of gluons from the physical part of the state space 
of QCD. For the scaling type of solutions with an infrared 
vanishing giuon propagator these violations can be shown 
analytically [7l|. This is not so for the decoupling type 
of solutions. Thus to investigate this question one has to 
determine the Schwinger function 

A(i) = Jd'xJ ^eM^P-^)D{p'), (48) 

numerically. Here D{p^) = Z(jj^)/p^ is the gluon prop- 
agator function. Our results for both type of solutions 
are shown in Fig. 1101 Apart from some variations in the 
scale this type of behavior can be seen for all our decou- 
pling solutions, i.e. it is independent of the value of G(0). 
This statement, however, may depend on the truncation 
scheme. 

From Fig. [To] it is plainly visible that, despite the dif- 
ferences in momentum space, the Schwinger function for 
both solutions is very similar. In particular, in both cases 
positivity is violated, and gluons are confined. This oc- 
curs for both cases at about the same scale of 1 fm, typical 
for the size of bound states. 

Furthermore, both results for the Schwinger function 
are in qualitative agreement with corresponding results 



from lattice calculations [9l|,|92]. Hence, despite the com- 
pletely different status with respect to the Kugo-Ojima 
framework, both solutions do not describe propagating 
gluons. 

In particular, the giuon is not characterized by a 
pole mass and an exponential decrease of a (positive) 
Schwinger function at large times. This implies that the 
mass parameter defined in (|42|) is also not a pole mass 
in the ordinary sense but at best a screening mass. In 
addition, a gauge-independent mass in a sense defined in 
(93! cannot be constructed. Hence, irrespective of its in- 
frared properties, the gluon is never an ordinary massive 
particle. 

IX. SUMMARY AND DISCUSSION 

We close with summarizing again some of our main 
results and conclusions. More detailed discussions have 
been provided in the related chapters. 

We have obtained, in various truncations of DSEs and 
FRGs, a one-parameter family of solutions for the ghost 
and gluon dressing functions of Landau gauge Yang-Mills 
theory. We argued that Slavnov- Taylor identities cannot 
be used to discriminate between these, contrary to what 
has been claimed in [20j . In particular, transversality is 
not indicative in this respect. Our numerical solutions for 
both, scaling and decoupling, are transverse on the level 
of numerical accuracy. However, global symmetries are 
indicative: exactly one member of this family, the only 
one exhibiting scaling behavior, is consistent with the 
existence of a globally well-defined BRST charge. The 
remaining solutions are of a decoupling type and are not 
BRST-symmetric. 

The appearance of the above mentioned one-parameter 
family of solutions in the DSEs and FRGs can be achieved 
by implementing boundary conditions 0, [lO] . The scal- 
ing solution is in accordance with the Kugo-Ojima and 
Gribov-Zwanziger confinement scenarios and solely gen- 
erated by the gauge fixing sector of the theory. In case 
of the decoupling solution the infrared dynamics is qual- 
itatively different. We have argued that the latter type 
of dynamics cannot be color-confining and preserving 
BRST symmetry at the same time. Both solutions, how- 
ever, satisfy the quark confinement criterion put forward 
in |33j . 

Finally, we would like to emphasize once more that 
neither solution corresponds to an ordinary massive gluon 
with a pole mass. Instead, in both cases the Schwinger 
function shows the presence of positivity violation and 
thus the gluon is not a free particle. 

In conclusion, each member of the one-parameter 
family of solutions is confining. 
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APPENDIX A: EXPLICIT EXPRESSIONS FOR 
THE DRESSING OF THE THREE-GLUON 
VERTEX AND THE INTEGRAL KERNELS 

The ansatz for the three-gluon vertex used in our trun- 
cation includes the dressing function 

with y = , z = (p — q)^ , the anomalous dimension 
6 = —9/44 of the ghost dressing function and the renor- 



malization constant Zi for the three-gluon vertex. The 
value of a is a free parameter which regulates the inter- 
action strength of the vertex. In Q a = 3(5 has been 
chosen. The scale Adec ~ 500 MeV is irrelevant for the 
scaling type of solutions since the gluon-loop including 
the three-gluon vertex becomes sub-leading at such small 
scales compared to the ghost loop. For decoupling type 
of solutions this scale implements decoupling within the 
three-gluon vertex. All results in this paper are qual- 
itatively independent of A^ec- Quantitative changes in 
the decoupling solutions are at the percent level when 
varying this scale within reasonable bounds. 



The kernels K, M and Q for the integral equations 
and p9p are given by 
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with the original kernels K, M and Q of the truncation scheme of ref. [sl. F70|: 
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APPENDIX B: MULTIPLICATIVE 
RENORMALIZABILITY OF THE GHOST DSE 

Multiplicative renormalizability of the ghost and gluon 
dressing functions implies the relations 

G(p^M')Z3(M^A2) = G°b^A2) 

z{p^^?)z,{^?A') - z\p\K^) 
Tf^\p^q^^,^)z^\^?A') = vf^^\p^q^t.^) (bi) 

between the bare dressing functions depending on a UV- 
cutoff A and denoted with superscript zero and the renor- 
malized dressing functions depending on the renormaliza- 
tion point p^ . 

Note that the right hand sides of eqs. (|Bip are inde- 
pendent of the renormalization point. Thus a finite re- 
normalization, i. e., a change in the renormalization point 



J 



from /i^ to is described by 



z{p'y) 



2 2^^3(A^^A2) 



Z{p\ii^) 



^3(/i^A2) 



Zg(y\k^) 



Furthermore we need the identity 



Z\ — ZgZy^ Zj,. 



(B2) 



(B3) 



Now let us consider the Dyson-Schwinger equation for 
the ghost propagator 



G{p\p?) 



Z3(M^A^)-Zl(/.^A2)g2(M^) 



7 (2^ q^ P ' ^ (P: 9, ) (B4) 



where P^^{k) 



Say 



denotes the transverse ant under a change of the renormalization point, since 



projector and = q^^ — p^^. This equation is invari- 



we have 
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1 ^ 1 ^3(^^A^) 



^3(A^^A2) 



d^q G{q\u^)Z,{v\K^)Z{e,U^)Z:,{v\k^) pT,,^p(2.1). 2^ ^l(M^A^) 
(2^)4 g2^3(A*^A2) A:2Z3(/i2,A2) ^A^^M^Wi. ^^''^'^ ^^^(^2^A2) 

A 

I 2 a2\ ^ / 2 \2\ 2/ 2\ f d 1 ) ^{^ ) r,T n \-n(2.1) / 2\ 

{v ,A )-Zi{v ,K )g {v ) ——J -2 r2 It^PuA^Wi 'iP,q,'^ ) 



(2^) 

ZJ{^,■^^')Z^gi,.^,A^)ZJ{l.^A^)Z■,iu\A^) 

Z2(i.2,A2)Z2(;,^A2)Z|(/.2,A2)Z3(M2,A2) 

A 

V/2a2\ ;?^2a2\2^2n/" "^9 ^(^ ' ^ ) ^(^ ) dT /; m^(2.1)/ 2n 

Z3(i^ ,A j-Zili' ,A)g {v) j J^lt^ -2 p Pt,Ak)Tl >{p,q,v ) 



(B5) 



(B6) 



(B7) 
(B8) 



This demonstrates explicitly that a solution G{p^ , i-i^) of 
eq. (|B4p and a solution G{p^ .v^) of eq. (|B8p are uniquely 

related by a multiplicative factor \ ■ Thus a 

Z3(l/-,A-) 

change of corresponding to g{^i^) g{v^) with 
/i^ must not generate a qualitatively different solu- 
tion. 

Since our truncation respects multiplicative rcnormal- 



izability this requirement is certainly respected. This is 
explicitly demonstrated in Fig. Illi where we compare 
our results for a{iJ?') = 1 and ^(jU^) = 0.5. The resulting 
running coupling is clearly an RG-invariant, whereas the 
ghost and gluon dressing functions arc related by finite 
re- normalization factors according to eq. (|B2p . as they 
should. 
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